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ABSTRACT 

Supersymmetric QCD corrections to top quark pair production by 77 fusion are 
calculated in the minimal supersymmetric standard model taking into account the 
effects of stops in the corrections to the total cross-section of tt production at the 
future e"*'e~ linear collider. We find that the relative correction can be a few percent 
for reasonable values of the parameters. 
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1. Introduction 

Recently, the top quark was found experimentally by the CDF and DO Collab- 
orations at Fermilab . The measured top-quark mass, 176 ± lOtll GeV, is close 
to the central value of that obtained from the best fit of the standard model (SM) 
to the latest LEP data. This is a remarkable success of the SM. However, there are 
a number of unsolved theoretical puzzles in the SM, and the latest LEP data on the 
branching ratio Rb oi Z ^ bb deviates from the SM prediction by 3.7 standard 
deviation[§]. These lead to more interest in considering possible new physics beyond 
the SM. Processes with top quarks may be good for testing new physics since the top 
quark is the heaviest particle yet found. Among various models of new physics so far 
considered, supersymmetry (SUSY) is a promising one at present. The simplest and 
interesting SUSY model is the minimal supersymmetric extension of the standard 
model (MSSM) 0. For solving the gauge hierarchy problem, SUSY should be bro- 
ken at energies around 1 TeV, and thus SUSY particles in the MSSM may be within 
the reach of future colliders. At the future e~^e~ linear collider with center-of-mass 
energy 500 GeV to 1.5 TeV, the e~^e~ —>■ tt event rate would be around 10^/yr, 
comparable with the Tevatron, however, the events would be easier to extract. It is 
possible to separately measure all of the various production and decay form factors 
of the top quark at the level of a few percent Thus theoretical calculations of the 
radiative corrections to the production and decay of the top quark is of importance. 
SUSY corrections to ti pair production in e~^e~ annihilation has been calculated in 
Ref.p]. At the e"*"e~ linear collider, hard photons can be obtained by laser backscat- 
tering. The intense 7 beams are generated by backward Compton scattering of soft 
photons from a laser of a few eV energy ||^. The luminosity distribution over the 

1 



77 invariant mass is broad and contains an abundant number of very energric pho- 
tons. The hard photon beam has approximately the same luminosity as the original 
electron beam. Therefore photon collisions are also good processes for testing new 
physics. 

In this paper we investigate the SUSY QCD correction to the top quark pro- 
duction by the process 77 ti in MSSM model. In Sec. II, we give our calculation 
of the SUSY QCD corrections to the scattering cross-section. The numerical results 
of the cross-section are given in Sec III. In recent years there have been renewed 
interest in the possibility of very light gluinos, with mass rrig < 5 GeV [0]. Our 
conclusion is that, with such light gluinos, the corrections can be large enough to 
be experimentally testable. 

2. The Cross-Section 

The total cross section of the production of ti in 77 collisions at the e'^e~ 
collider is obtained by folding the elementary cross-section for the processes 77 — > tt 
with the photon luminosity ||^ 



as = I dz — r^o"(77 tt at s = z'^s), (1) 

J2mt/\^ dz 

where ^/s{V§) is the e"'"e~(77) center-of-mass energy and the quantity ^^j^ is the 
photon luminosity defined as 

^ — F,/e(x)F^/,(zVx). (2) 



dz 

For unpolarized initial electrons and laser,the energy spectrum of the back-scattered 
photon is given by 

1 1 Ax Ax"^ 

File = :d(^[i-^ + 1^-^(13^ + ^(1^^^ 



where ^ = 4EoUJo/ml with rrie and Eq the incident electron mass and energy, re- 
spectively, and the laser-photon energy, x is the fraction of energy of the incident 
electron carried by the back-scattered photon. Following Ref.[||, we choose ^ and 

•^max be 

e = 2(1 + V2) ^ 4.8, x^,, ^0.83, D(0 ~ 1-8. (5) 

The Feynman diagrams contributing to the SUSY 0{asOte) corrections are shown in 
Fig.l. In our calculation, we use dimensional regularization and take the on-shell 
renormalization scheme. 

In Fig.l one has to include the contributions of both stops. As is well-known 
[P], the supersymmetric partner of left- and right-handed massive quarks mix with 
each other. The mass eigenstates qi and q2 are related to the current eigenstates ql 
and qR by 

qi = qiCOsOq + qusinOq, q2 = -qisinOq + qRCOsOg. (6) 

The mixing angle 6q as well as the masses m^-^ , m^^ of the physical stops can be 
calculated from the following mass matrices [1^ 

/ 



^ —mt{At + ficotp) mf^ + mf + 0A6Dz 



(7) 



where Dz = MzCos2f3, tan(3 is the ratio of the vacuum expectation values of the 
two neutral Higgs fields of the MSSM, m^-^, m^-^ are soft breaking masses. At are pa- 
rameters describing the strength of nonsupersymmetric trilinear scalar interactions, 
and /i is the supersymmetric Higgs(ino) mass which also appears in the trilinear 
scalar vertices. 



In the presence of squark mixing, the squark-quark-gluino interaction La- 
grangian is given by 

Lgqq = -iV^gsT''q[{cos9gqi - sin9gq2)^^ 
-{sin9gqi + cos9gq2)^^]ga + h.c, 

where gg is the strong couphng constant and are SU{3)c generators. 

To the precision of the ©(a^oje) SUSY corrections, the renormahzed amphtude 
for 77 —>■ tt is 

Mren = Mq + 5M'^^^ + SM"^''^^'' + 5M^°'' + 5M^ (9) 

where Mq is the tree-level amplitude and 5M represents SUSY QCD correctrions, 
which are 

Mo = e>'{PA)e''{p^)u{p2),T^MP^) W 



with 



and 



T^v = 7^^57,.(;^3 - A + mt)7^ + --^^^7^(^4 - A + mt)7^, (11) 

b //t-^ Uj //tf 



5M'^^^ = 5M^^'-^(*) + (5M*"'-^(") (12) 

^ -^vertex ^ ^j^vertex{t) ^j^^vertex{u) ^g-j 
5M^°^ = 5M''°^(*) + 5M'"'^("\ (14) 



(5M-'/(*) = e/'(p,)e-(p3)^z(p2)^^^7^(^2-A + ^t)(-«^) 
x(A - A +"^t)7i^^^(Pi), 

x(-iQ?7.) + (-^Qf7^)(A - A + ^t)^Af )]}^(pi). 



(15) 



(16) 
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, nbox(t) , nbox(t) , nbox(t) , I nbox(t) 

+PltiP2uIs + P2t,Pluf9 + P2,,P2uflO + Piluli^Ju 

, I rbox(t) , I rbox(t) , I rbox(t) 

, I rbox(t) , I rbox(t) , I rbox(t) 

+P4P2ulf^fl5 +P4P2f,lufw +PiPlpiPlvJn 

+hPii^P2vf\T^^ +i)iP2tiPivf\T^^^ + i>mi^P2vfll''^^}v{pi), 

SM^ = tMt^{p;)u{p2)QlSv{p,), 



(17) 



with 

-it = /f^n//?*^ (19) 

" a'"' a'"' a(") 

^Af) = 7m/i" +P2m/2^ +M2,ft (20) 

" /_ir,N A ('2) A ('2) A ('2) 

^Af^ = luft-' +Piuft +hPiuft , (21) 

where Qt = 2/3, S and the form-factors fi are given in the Appendix. Instead of 
calculating the square of the amplitudes explicitly, we calculate the amplitudes nu- 



merically by using the method of Ref . |]TT| . This greatly simplifies our calculations. 



We only explicitly give the results for the t-channel contributions to the SUSY 
corrections. The u-channel results can be obtained by the following substitutions 

P3 ^ P4, ^ T\ i ^ u. (22) 



3. Nemerical Results and Conclusion 



Now we present the numerical results. We take rrit 
two-loop running coupling constant and the input cte 
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= 176 GeV, and use the 
= 1/128. For the SUSY 



parameters involved in our calculations, we see that once tan(3 and nif^ are fixed, we 
are, in general, free to choose two independent parameters in the stop mass matrix, 
namely nif^ and At + ficotp, or equivlantly nif^ and m^^. To avoid the singularities 
at small angles, we take the following kinematical cuts 

\ri\ < 2.5, pt > 20 GeV. (23) 

In this kinematical region the relative corrections are actually large. In the numerical 
calculation, we have checked our program with the requirement of gauge invariance 
to the accuracy 10^^°. 

In Figs. 2-8 we give the numerical results in a simple case in which we set 
tanP = 1 and rrif^ = rrif^ = rriq (coresponding to the mixing angle equal to vr/4). 
We choose m^- as the light stop mass and require it to be heavier than 45 GeV . 

Fig.2 (Fig. 3) show the dependence of the corrections on nif^ for fixed nig = 
3 GeV, nig = 150 GeV (450 GeV) and = 0.5 TeV (1.5 TeV). We see that the 
corrections can be either positive or negative depending on the light stop mass. The 
corrections become their negative maxima at m^- = 170 GeV. Fig. 4 (Fig. 5) show the 
dependence of the corrections on for fixed rrig = 3 GeV, rUg = 450 GeV and m^^ = 
50 GeV (150 GeV). The corrections are positive when m^^ = 50 GeV, and negative 
when = 150 GeV. When ^/s varies from 0.5 TeV to 1.5 TeV, the corrections 
vary from 2% (-3.6%) to 0.9% (-3.8%) in the case of m^-^ = 50 GeV (150 GeV). 
Fig. 6 (Fig. 7) shows the dependence of the corrections on rriq for fixed rrig = 3 GeV, 
m^^ = 50 GeV (150 GeV) and = 0.5 TeV (1.5 TeV). The relative corrections 
increase (decrease) as rriq varies from 50 GeV to 450 GeV in the case of m^^ = 
50 GeV (150 GeV). The largest relative correction in Fig. 7 can exceed —5%. Fig. 8 
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shows the dependence of the corrections on rrig for fixed m^^ = 150 GeV, rriq = 
450 GeV and ^/s = 0.5 TeV (1.5 TeV). If the glunios are hght§, e.g. m-g = 3 GeV, 
the corrections can reach —4%. Whereas for heavy glunios as rrig > 100 GeV, the 
corrections are less than 1%. 

We have also done the numerical calculations for tanf3 = 10 and found that the 
corrections are not sensitive to the value of tan(3. For example, with nig = 3 GeV, 
rui^ = 150 GeV, nig = 450 GeV and y/s = 0.5 TeV, we get Aa = -3.60 for tan(3 = 1 
and Act = —3.58 for tanf] = 10. 

We conclude that, in the case with stop mixing and light gluinos like nig 
around a few GeV, the SUSY QCD corrections to the cross-section of top quark 
pair production in 77 fusion at the e^e~ collider can be as large as a few percent of 
the tree-level cross-section. This is experimentally testable. 
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Appendix 



We give here the form factors for the matrix element. They are written in 



terms of the usual one-, two-, three- and four -point scalar loop integrals of Ref . |[T2 1 



ai = -62 



fl 



—j={cos9 — sin9), 02 = &i = -y={cos9 + sini 



Arv cv ■ 

Y.[iF^B, - m,ZV, - DM,)] 



3 

Aasae 



J2[G,{Bo + B,) + ZVi], 



where Bq, Bi{t,m£^,mg) are 2-point Feynman integrals |T2|. 



—i- 



f7 



x(tl) 



3 

I 

3 

.4aoap 



i 

^[2G,mi(Cn + C21) - 2F,iCo - Cu)] 

i 

J2[2G,{-C,2-C23)], 



where Go,Gij{—p2,P4,,mg,'m£.,m{J are 3-point Feynman integrals 

.a(*2) AasQe 



fl 



h 



i(t2) 



— l- 



3 

AOLsOie 
AdsCle 



Y.[2G,G2^ + ZV,] 

i 

Y,[2Gmt{-Cu - C2l) + 2F,(Co + di)] 

i 

^[-2G,(Ci2 + C23)], 



where Go,Gij{pi, —ps,rng,m£-,m{.) are the 3-point Feynman integrals 



/r 

/2 



box{t) 



^CisOte 
Aci^df 



Y,[2FiD27-2GmtD:,ii] 

i 

J2[2FiD2r-2G,mtDsn] 



(24) 
(25) 

(26) 
(27) 



(28) 
(29) 
(30) 



(31) 
(32) 
(33) 



(34) 
(35) 



box{t) 



box(i) 



fbox{t) 



fbox{t) 
J6 
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box{t) 



— I- 



fi 



box{t) 



^ i 

AOLgOLe 
I 

3 

AoisOe 



— I- 



E[4<^«^313] 

i 

5:[4G',(-Z^3ii+^3i2)] 

i 

^[4G',(-L>27-i^311+i^313)] 

i 

^[4Fi(L>i3 + L>26) - AGimt{D^^ + L'gio)] 

i 

^[4F,(-L>25 + D^^) + 4G',mt(L>35 - L'aio)], 



box{t) 



•4asae 



E[4Fi(-L>o - ^11 + ^13 - - D24 + + 4G',x 



nbox(t) 

JlO 



E[4F,(L>n - + D21 - D24 - D25 + D2e) + AGiX 



mt{-D2i + D24 - D31 + D34 + L'35 - ^310)], 



nbox(t) 
111 



rbox{t) 
JI2 



Y,[^Gi{D^^2 - i?313)] 

i 

Y,[^G,{D^^2 - i^313)] 



nbox(t) 

/l3 



= 



fbox{t) 
Jl4 

nbox{t) 

Jl5 








rbox{t) 
/16 

pbox{t) 
111 



pbox{t) 
/18 







4q; Oi 

^^r^ E[4Gi(-^23 + D2& + L»38 - L'39)] 



3 

4Q;<,Q;f, 



= I- 



310) 



f. 



box{t) 
19 



E[4G'i(L'i3 - i?12 - ^23 - ^24 + D. 



25 



+2^26 - ^22 - i^36 + ^38 - ^39 + ^31o)], 



(53) 

-D35 - + 1^37 + ^38 - ^39)], 



S = i-asaeY^[GimtCu- FiCo]{-p2,P4+P3,rng,mf^,mfJ, (54) 



where Dq, Dij, Dijk{~P2,P4,P3,'rng,m£.,m£^,m£.) are 4-point Feynman integrals ||T2 

The renormahzation constants are 

ZV, = -G,[Bo + B,]{p,m£,mg)\p2^^2-[2m^,G,-^{Bo + B,) 



9p2. . . ^^^^ 

T-n ~ TOO ~ I „ n 

DM, = [FiBo + mtG,{Bo + Bi)]{p,m£,mg)\p2^^2 (56) 



-2mtFi^Bo] (p, m£ , m-g) \p2._ 
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Figure Captions 

Fig.l Tree-level Feynaman diagrams and 0{asae) SUSY QCD corrections to 77 — > 
tl 

Fig. 2 Relative SUSY QCD corrections to 77 fusion cross-section versus rrif^ for 
mg = 3 GeV, niq = 150 GeV. 

Fig.3 Same as fig.2,but for = 3 GeV, mq = 450 GeV. 

Fig.4 Same as Fig. 2, but versus ^/s for = 3 GeV, mj^ = 50 GeV, mq = 450 GeV. 
Fig.5 Same as Fig. 2, but versus for nig — 3 GeV, mj^ = 150 GeV, mq = 450 GeV. 
Fig.6 Same as Fig.2, but versus rUq for rrig — 3 GeV, mj^ = 50 GeV. 
Fig.7 Same as Fig.2, but versus rriq for rrig — 3 GeV, mj^ = 150 GeV. 
Fig.8 Same as Fig.2, but versus nig for niq — 450 GeV, m^^ = 150 GeV. 
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